Universality and itinerant ferromagnetism in rotating strongly interacting Fermi gases 
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We analytically determine the properties of three interacting fermions in a harmonic trap subject to an external 
rotation. Thermodynamic quantities such as the entropy and energy are calculated from the third order quantum 
virial expansion. By parameterizing the solutions in the rotating frame we find that the energy and entropy are 
universal for all rotations in the strongly interacting regime. Additionally, we find that rotation suppresses the 
onset of itinerant ferromagnetism in strongly interacting repulsive three-body systems. 
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Broad Feshbach resonances in two-component atomic 
Fermi gases have made it possible to explore the crossover 
between Bardeen-Cooper-Schrieffer (BCS) superfluidity and 
Bose-Einstein condensation (BEC) |1-3 |. In the strongly in- 
teracting regime the difficulties in developing many-body the- 
ories for these systems has motivated the study of exact solu- 
tions to few-body problems as a means to gain insight into the 
many-body problem B HTTI . Of particular interest has been 
the prediction [[TOl fT2UT4l and observation [fTSl 161 of uni- 
versality in the strongly interacting regime of thermodynamic 
quantities, such as the energy and entropy. To date such cal- 
culations have been restricted to gases in non-rotating spheri- 
cally symmetric harmonic traps. 

In this work we consider the propeities of a strongly inter- 
acting gas in a rotating trap. In particular we solve the two- 
IfTTl and three-body problems and calculate the virial expan- 
sion of the thermodynamic potential to third order, enabling 
the calculation of thermodynamic quantities. Interestingly, we 
find in the thermodynamic limit that the second- and third- 
order virial coefficients are universal, with respect to an exter- 
nal rotation and trapping frequencies. From this we show that 
thermodynamic quantities such as the energy and entropy are 
universal with respect to rotation through a simple rescaling 
of the Fermi energy. 

In addition to these thermodynamic results we also exam- 
ine the interplay between rotation and the emergence of itin- 
erant ferromagnetism in strongly interacting ultra-cold Fermi 
gases. In the original work of Stoner |18| it was proposed 
via a meanfield theory that a repulsive Fermi gas will al- 
ways exhibit a ferromagnetic phase. Most recent experimen- 
tal evidence |19| suggests there is no transition. Current 
theoretical work, Monte Carlo simulations and Tan relations 
Hoi I20h22|| are contradictory. In this work we show that for 
the three-body problem, rotation suppresses the emergence of 
itinerant feiTomagnetism. Additionally, in the thermodynamic 
limit itinerant ferromagnetism is suppressed for temperatures 
T > 10^^ Tp as the rotation frequency approaches the trapping 
frequency. 

Our starting point is the wavefunction, i/'(ri,r2, . . .), of 
particles interacting in the s-wave channel at low energies. 
This satisfies the Bethe-Peierls boundary conditions 



lim dirijtf/)/dnj = -nji///a, 

n,~»0 



(1) 



a, and = |r,-rj| is the separation of opposite-spin fermions. 
Away from r,j - 0, the wavefunction of A^ particles in a spher- 
ically symmetric rotating harmonic trap satisfies the non- 
interacting Schrodinger equation. 
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where r, and jj. are the position and mass of each particle and 
o) and are the trapping and rotation frequencies, the latter 
assumed to be defined about the z axis. 

The center-of-mass Hamiltonian can be decoupled from 
Eq. ^ and defines the rotating harmonic motion of a parti- 
cle of mass M - Nfi with energy spectrum 



(2n + l + 3/2)na) + mtiQ.^ 



(3) 



where n, I, m label the usual harmonic oscillator eigenstates 
R„iY'". The relative energy, £^1 - E - E^^^, incorporates 
the effects of the contact interaction but not the effects of the 
external rotation, since only s-wave states may interact. For 
two opposite-spin fermions the wavefunction in relative co- 
ordinates that satisfies the Bethe-Peierls boundary condition 
Eq. ([TJ is a 

i/r5,'(r; v) oc r (-v) U (-V, 3/2, r^jcf) exp [-r^lcF] , (4) 

where U is the confluent hypergeometric function of the sec- 
ond kind. A pseudo-quantum number, v, parametrizes the rel- 
ative energy E^^x - (2v + 3/2)^w, and satisfies the relation 



2r(-v) 
r(-v-i/2) 



(5) 



where the interaction is parametrized by the scattering length 



for harmonic oscillator length d - ^jtiJiJujS). In particular, for 
the unitary limit, where a +oo, the relative energy spectrum 
simplifies to Ei-ei = (2n + 1 /2)fujj, where n is any no-negative 
integer. 

To find the energy spectrum of three interacting fermions 
in a rotating trap, we consider the configuration of two 
spin up fermions and one spin down, tit- where two 
opposite spin particles interact at a point and form a 
pair, and the third moves relative to the pair. We de- 
fine the center-of-mass coordinate of the three paiticles as 
R = (ri H- r2 H- r3)/3, the relative coordinate between the in- 
teracting pair, r = ri - and the relative coordinate between 
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the third non-interacting particle and the center-of-mass of the 
pair as p = (2/ V3)[r3 - (ri + r2)/2]. In this Jacobi coordi- 
nate system the center-of-mass Hamihonian decouples from 
the relative Hamiltonian 



(6) 



where /i/2 is the reduced mass of the interacting pair. Like 
the two-body system, the angular momentum vanishes in the 
interacting pair but the third fermion can rotate around the 
pair and be affected by the external rotation, Q^. This couples 
higher order angular momentum states to lower energies in the 
system. In order to solve the relative Hamiltonian (|6]l we take 

oo 

'V^l{V,p) = (1 - Pn) 2 C„<AS!(r; VnUn)Rnm{p)YTip\ (7) 

as an ansatz for the wavefunction, where f 13 is an operator 
that exchanges the spin t particles. The eigenenergies of this 
system are 



El-el - 



tia> + mfiQ... 



(8) 



The presence of the rotational term, mMl^, in the eigenen- 
ergy spectrum shifts the non-rotating energy spectrum found 
in f9 |. To solve for the coefficients c„ in Eq. (|7]) we use the 
Bethe-Peierls boundary condition Eq. ([TJ and choose a set of 
quantum numbers nlm and energy E^^\ to solve for a particular 
v„i,„ and scattering length a Q- 

The relative energy spectrum can be found numerically for 
any scattering length, but in the unitary regime it is simpler 
to use the method of Werner and Castin IS) to obtain the 
energies using hyperspherical coordinates (R,a,r,p), where 
R - yjir^ + p^)/2 is the hyperradius, a - arctan(r/p) is the 
first hyperangle and r and p are the direction of each Jacobi 
coordinate. Using this coordinate system and the ansatz for 
the wavefunction |6], 



sin(2a') 



(9) 



the Hamiltonian, Eq. (|6]), can be written as two decoupled 
Schrodinger equations. 



\F" + -F'\ + 

2m\ R 



^ H — mo) R — mTiD., 
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and 



/(/ + 1) 2 
cos^(a) 



(10) 
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For three fermions is always positive and we can interpret 
the hyperradial Schodinger equation ( [TO] l as a particle mov- 
ing in a two dimensional effective potential {tfs^i/(2mR^) + 
^maP'R^) with energy spectrum Pl 



Eiei — {2q + Sni + lytui) + mhO,^, 



(12) 




FIG. 1. Energy spectrum of three interacting fermions with / = 2 and 
m = -2, and 2 for plots (a), (b) and (c) respectively with a rotation 
of flj = 0.9a). We can see the shifting of the energy spectrum with 
the rotation included, in particular we can see the number of lower 
energy states increase for m = -2. 



where ^ is a positive integer The solutions to Eq. ([TT} must 
satisfy (/?(;r/2) = so that the ansatz Eq. (|9]) does not diverge. 
Hence, 



(p(a) oc cos(a) 



/+lp(/+l/2,-l/2)r 

(.v-/-i)/2 L 



- cos(a)]. 



(13) 



where P'"J'^\x) is the regular Jacobi polynomial [23 1. The 
eigenvalues are determined from the Bethe-Peierls bound- 
ary condition ([T]), which in hyperspherical coordinates reads 



,,„,_(_,yJ_,(f) = „ 



(14) 



The procedure for solving Eq. ( [T4] i for values of s„i using the 
general solutions to the hyperangle equation is given in [10|. 
It can be shown that the two spectra Eq. ([12]) and (|8]l are the 
same in the unitary limit, where the solutions from Eq. (jSj are 
found numerically. 

The energy spectrum of three interacting fermions is shown 
in Fig. ([1} for a rotation of Q- = 0.9a>, a relative angular mo- 
mentum of / = 2 and axial angular momentum of m = -2, 0, 2. 
We see the energy levels for the three interacting fermions 
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FIG. 2. (a,b), Energy per particle E/(NEp) and (c,d), entropy per par- 
ticle S /(NkB) for dimensionless rotation ^ as a function of reduced 
temperature T/T^, in the strongly (a,c), attractive and (b,d), repul- 
sive regimes. For comparison we plot the temperature at which the 
second and third order virial expansions differ by 1% for a rotation 
^, the dashed line. The hashed region in plots (a) and (c) indicate an 
unphysical solution to the virial expansion. 



shift for each m quantum number, lifting the (21 + 1) fold de- 
generacy. 

To obtain the repulsive spectrum we omit the solution en- 
ergy levels which are the lowest in each /-subspace lowest or- 
der n - bound state in each / subspace, i.e. the sqi energy 
levels. The lowest energy in the repulsive regime is then the 
relative energy, which is dependent upon the rotation Q.^, plus 
the center-of-mass energy, 

£]f ^ (^^^^ + + ,„;^Q_ + l_5nco. (15) 

This can be compared to the energy of three polarized 
fermions, which is the sum of the three lowest allowed en- 
ergies: 

£JP = i.Shoj + i.Shoj + [i.Shoj + n,)n]. (16) 

Comparing Eqs. ( [T5| l and ([T6| it is possible to find a critical 
rotation, Qj., for which the ground state energy of three re- 
pulsively interacting fermions becomes lower than three non- 
interacting polarized fermions. This is a rotation for which 
the three non-interacting fermions are unstable with respect 
to the three interacting fermion system. Choosing the lowest 
relative energy level in the / - 2 subspace, si2 - 4.80, and a 
magnetic quantum number m = -I, we find a critical rotation 
of D.C - 0.8w. Since the ground state energy of three repul- 
sively interacting fermions can be controlled by varying an 
external rotation, this indicates that an itinerant ferromagnetic 
transition, in the three-body system, only occurs for a rotation 

a < Qc. 

The two- fl7] and three-body solutions can now be used 
to calculate the many-body properties of strongly interacting 
rotating Fermi gases. This is done via a quantum virial expan- 
sion of the grand thermodynamic potential, O - -ksT In X, 



in terms of the fugacity z: 

1 r" 

2 Jo 

-kBTQi{z + Ab2Z^ + Ahz' + ...), (17) 

where 

Ab2^AQ2/Qu (18) 
Abi^AQ,/Qi-AQ2, (19) 

and AQ„ = Q„ - with Qn = Tr [exp(-'7/^/A:Br)] r24|. To 
first order for dimensionless rotation ^ = Q./w, Qi is given 
by 

Following IfTTll . the second-order virial coefficient for a rotat- 
ing trapped gas in the high temperature hmit is 

Atf ^\-'^ + ... (21) 
2 4 32 

A^- = -i + ^ + ... (22) 

where a> - Hot /k^T is the reduced trapping frequency and 
(i) — » represents the thermodynamic limit. Extending the 
work of non-rotating systems in |9 | we find that Ab"^" is uni- 
versal for any rotation and given by 

Abf ^ -0.06833960 + 0{o?). (23) 

For a repulsive gas the third-order virial coefficient to lowest 
order in Co is also universal, 

A/77 ^ 0.34976 -I- 0{(h). (24) 

Despite the rotational dependence of the two- and three-body 
eigenspectmms the second- and third-order virial coefficients 
are independent of rotation in the thermodynamic limit. 

We are now able to calculate the total energy E - -3<E) and 
entropy S - -d^/dT from the thermodynamic potential O of 
a strongly interacting gas [24 1. To determine the thermody- 
namic potential, the fugacity z = expQ^ /k^T) must be calcu- 
lated from the total number of particles = -cJO/S//, where 
/i is the chemical potential. 

In Fig. I2] we plot the energy, in units of NE'^, where 
Ep — QNy^Tiu is the non-rotating Fermi energy, and entropy 
per particle as functions of reduced temperature, - E^jk^, 
and rotation using the virial expansion to third-order for a 
strongly repulsive, (b,d) and attractive, (a,c) Fermi gas. The 
hashed areas in plots (a) and (c) correspond to solutions of the 
energy and entropy that are unphysical. The dashed curve in 
Fig.|2]is the temperature at which the second and third order 
expansion differ by 1%, providing a conservative estimate of 
the temperature range of validity for the virial expansion. As 
the rotation is changed, for a fixed temperature, the energy 
(a,b) and entropy (c,d) change. Hence, the thermodynamic 
quantities appear not to be universal with respect to rotation. 
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FIG. 3. Chemical potential, (a), entropy per particle, (b) and energy 
per particle, (c) of a strongly attractive (solid), repulsive (dashed) and 
ideal Fermi gas (dotted) for dimensionless rotation f = 0, 0.5, 0.9 and 
0.99 as a function of reduced temperature T/rf . 



From the three-body calculations Eqs. ( [T5] l and ( fTS] ) we 
see that itinerant ferromagnetism is suppressed for Q- > O.Sw. 
Fig. |2|b) plots the total energy in the strongly repulsive 
regime. As ^ — » 1 the validity of the solutions extends to 
r — > 0. In this regime, by comparing the energy of the 
strongly interacting gas with the equivalent non-interacting 
polarized gas, we find that the itinerant ferromagnetic phase 
is suppressed for r/r" = 10"^ |25|. 



Figure |2] demonstrates that for a given rotation frequency, 
^, there is a universal dependence of E/(NE^) and S /(Nks) 
as a function of Tp, with respect to particle number and trap- 
ping frequency. However, through a simple rescaling of the 
Fermi energy [temperature] of the form - £"(1 -^^)"'^-' 
[Tp - Tp^l -^^)"'^-'] it is possible to remove the rotational 
dependence observed in Fig. [2] This generalized universality 
arises from the fact that under this rescaling the functional de- 
pendence of ^ is removed from both the single particle cluster 
function, Qi [Eq. (|20]i], and the chemical potential. To em- 
phasize the universal nature of the chemical potential with re- 
spect to rotation, fi/E^ is plotted in Fig. 3 [a) as a function of 
T/T^ for various rotations, ^, in the strongly attractive (solid), 
repulsive (dashed) and ideal (dotted) regimes. Hence, in con- 



junction with the fact that the second- [Eqs.(21|22i] and third 



order [Eqs. (23124 1] virial coefficients are independent of ro- 



tation in the thermodynamic limit, the thermodynamic poten- 
tial is independent of rotation. As a direct consequence the 
rescaled energy, E/iNE^), and entropy, S/iNks), are univer- 
sal with respect to rotation as functions of T/T^,. This prop- 
erty is demonstrated in Figs. |3jb,c) which plot E/(NE^) (b) 
and S /(NkB) (c) as a function of T/T^ for various values of ^ 
in the strongly attractive (solid), repulsive (dashed) and ideal 
(dotted) regimes. 

In conclusion we have examined the problem of three 
ultracold fermions in a harmonic trap subject to an external 
rotation. For the three-body problem we have demonstrated 
that rotation suppresses the transition to a ferromagnetic state. 
Explicitly we have shown that the three-body ferromagnetic 
state has a higher energy than the strongly interacting re- 
pulsive state for > 0.8a; and is, consequently, unstable. 
Additionally, from the three-body solutions and the use of pre- 
vious two-body results |T7l we have calculated the equations 
of state using the virial expansion to third-order. Despite the 
rotational dependence of the two- and three-body eigenspec- 
trums we have generalized the universal nature of strongly 
interacting fermions to include rotation by a simple rescaling 
of the Fermi energy and temperature. These results could be 
used as benchmarks in experiment to test the universal prop- 
erties of strongly interacting rotating ultracold Fermi gases. 
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